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Abstract 

We present clustering measurements for samples of galaxies selected by morphological type and 
luminosity from the recently completed Stromlo-APM Redshift Survey. We find very different results 
between real and redshift-space estimates of the correlation function. The real space correlation 
function for the all-galaxies sample is well fit on scales 0.2-20/i~"^Mpc by a power-law with slope 
7j. = 1.71 and correlation length tq = 5.1/i~-^Mpc. In redshift space the slope is shallower, 7^ = 1.47 
and the correlation length is slightly higher, sq = 5.9/i~"^Mpc. 

Early type galaxies are clustered more strongly by a factor 3.5-5.5, than late type galaxies. In real 
space the slope of the correlation function for early type galaxies is 7^. = 1.85, slightly steeper than for 
late types, 7^. = 1.64. In redshift space however, early type galaxies have a very shallow correlation 
function slope, 7^ = 1.25. This implies that these early-type galaxies suffer from enhanced redshift- 
space distortions compared to late-type galaxies. 

Low-luminosity galaxies are clustered more weakly by a factor of ~ 2 than L* and brighter galaxies 
on scales > l/i~-^Mpc. Also the slope of the correlation function is steeper for low-luminosity galax- 
ies, so that the amplitude is a factor 4 lower at 10/i~-^Mpc. No difference, however, is seen between 
the clustering of X* and more luminous galaxies, an observation which may be hard to reconcile with 
some theories of biased galaxy formation. Both redshift-space and real-space clustering estimates 
show a similar dependence on luminosity. 

Our results hint that luminosity segregation may be a purely primordial effect, due to a lower bias 
factor for low-luminosity galaxies, whereas morphological segregation, being most apparent on scales 
< l/i~-^Mpc, may be enhanced by environmental factors. 



Subject headings: galaxies: clustering — galaxies: fundamental parameters — large-scale structure 
of Universe 
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1 Introduction 



The dependence of galaxy clustering on luminosity and morphology is of great relevance to all galaxy 
formation theories and in understanding the large-scale matter distribution in the universe. If there 
are differences between the clustering of various types of galaxies we can immediately infer that at 
least one of the galaxy types is a biased tracer of the underlying mass distribution. In general we 
expect all galaxy samples to be biased at some level relative to the mass, and the differences in 
clustering can test various models for the bias of galaxies relative to mass. For example, the process 
of "natural bias" (White et al. 1987) leads to a galaxy correlation function that is a constant factor 
times the mass correlation function, with the factor being larger for more massive galaxies. More 
realistically, galaxy formation depends on complex processes which involve the local environment, 
some feedback mechanisms or galaxy interactions, as well as the depth of the dark matter potential 
{e.g. Dekel and Rees 1987). Quantitative measurements of the relative distribution of galaxies of 
different luminosities and types will tightly constrain models of these processes. 

That the correlation function for elliptical and lenticular galaxies has a steeper slope and larger 
amplitude on small scales than the correlation function for spiral galaxies has been known at least 
since Davis and Geller (1976) calculated the angular correlation function of morphologically-selected 
galaxy subsamples from the Uppsala Catalogue (Nilson 1973). The possible dependence of galaxy 
clustering with luminosity has been much more controversial. Several groups {e.g. Bothun et al. 
1986, Binggeli et al. 1990, Eder et al. 1989, Thuan et al. 1991, Weinberg et al. 1991) have found 
that the clustering of dwarf galaxies is consistent with that of bright galaxies (or at least that dwarf 
galaxies do not 'fill the voids' in the bright galaxy distribution). Alimi et al. (1988) and Phillipps and 
Shanks (1987) have measured clustering strength over a range of luminosities and find no evidence 
for enhanced clustering of more luminous galaxies. On the other hand, Davis et al. (1988), Hamilton 
(1988), Salzer et al. (1990), Santiago and da Costa (1990), lovino et al. (1993) and Park et al. 1994 
do claim to detect stronger clustering of bright galaxies compared with faint. Maurogordato and 
Lachieze-Rey (1991) find luminosity-dependence in the void probability function but not in the two- 
point correlation function. Hasegawa and Umemura (1993), after extinction-correcting CfA data, 
find weak luminosity segregation of opposite sign in early and late type galaxies. 

All of the above analyses use fairly shallow (miim ~ 15) catalogues of galaxies so that intrinsically 
faint galaxies can only be seen within a very small, nearby volume. It has thus been difficult to 
make a reliable comparison between clustering properties of bright and faint galaxies. In this paper 
we analyse the Stromlo-APM Redshift Survey which samples a much larger volume of space than 
any previous optically-selected galaxy redshift surveys; L* galaxies can be seen out to a distance 
~ 180/i~"^Mpc, that is within a volume ~ 2.5 million h~^Mpc^. This large volume has been surveyed 
rapidly by using a 1 in 20 sparse sampling strategy to select galaxies for redshift measurement. The 
uniform sampling of such a large volume makes the Stromlo-APM Survey an extremely powerful 
sample for studying luminosity segregation in galaxy clustering. 

The construction of the survey has been briefly described in an earlier paper (Loveday et al. 1992a, 
hereafter Paper 1) and will be described in full in a future paper in this series. An analysis of the 
large-scale clustering of galaxies in the Stromlo-APM Redshift Survey and a comparison with two 
versions of the Cold Dark Matter theory has been given by Loveday et al. 1992b (hereafter LEPM). 
In the present paper we study the clustering properties of subsamples of galaxies selected from the 
survey by morphological type and luminosity. The paper is organised as follows. In section 2 we 
describe the galaxy samples used in the clustering analyses. In section 3 we discuss two estimators 
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for the redshift-space correlation function ${s), one dependent on and one independent of an assumed 
mean galaxy density and apply these estimators to the galaxy samples listed in section 2. In section 
4 we present two estimates of the real-space correlation function ^(r) unaffected by redshift-space 
distortions. Our conclusions are presented in section 5. Throughout the paper, we use r to denote 
real-space separations and s to denote separations in redshift-space. Unless otherwise stated, error 
bars in figures and quoted errors in numerical quantities are 1-sigma dispersions calculated by 
analysing nine bootstrap-resampled versions of the survey (Barrow et al. 1984). 

2 The Galaxy Samples 

The Stromolo-APM Redshift Survey consists of 1787 galaxies with bj < 17.15 selected randomly 
at a rate of 1 in 20 from the APM Galaxy Survey (Maddox et al. 1990a, b). The survey covers 
a solid angle of 1.3 Sr (4300 square degrees) in the south galactic cap. The APM magnitudes have 
been calibrated and corrected for photographic saturation using CCD photometry as described in 
Paper 1. An approximate morphological type was assigned to each galaxy by visually inspecting 
the images on the UKST survey plates. Redshifts have been obtained with the MSSSO 2.3m 
telescope at Siding Spring. Measured radial velocities are transformed to the local group frame 
using ■;; = ■!; + 300 sin(Z) cos(6) and we assume A = 0, qo = 0.5 and Hq = 100 km s~-^Mpc~-^ with 
uniform Hubble flow in calculating distances and absolute magnitudes. We adopt A;-corrections for 
different morphological types in the bj system as described by Efstathiou, Ellis and Peterson (1988). 
More details about the survey are given in Paper 1, and the construction will be described in full 
in a future paper in this series. 

We draw six samples from the Stromlo-APM Redshift Survey: (a) all galaxies; (b) low, (c) medium 
and (d) high luminosity galaxies; (e) early and (f) late morphological types. These samples are 
defined in Table 1 and in Fig. 1 we show the redshift-distance histogram for each sample. Note that 
the volumes of each sub-sample all overlap to some extent, and that even the faintest sample extends 
beyond 100/i~-^Mpc in depth. For the luminosity-selected samples (b, c and d) we have applied 
an apparent magnitude bright limit, m > 15, since galaxies brighter than 15th mag suffer from 
photographic saturation on deep Schmidt plates and hence unreliable magnitudes. The luminosity 
limits were chosen to divide the galaxies into sub-X*, ~ L* and super-X* samples with roughly equal 
numbers per sample. 

3 The Redshift Space Correlation Function ^(s) 
3.1 Estimators 

In any flux-limited survey, the observed density of galaxies will decrease with distance x from 
the observer. In order that an estimate of the correlation function not be dominated by the nearby 
galaxies, it is important to give the appropriate weight to each galaxy. The variance in the estimated 
^(s) is minimised by weighting each galaxy in a pair at redshift-space separation s as 




(1) 
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(See Appendix; this weighting scheme was first used by Efstathiou (1988)), where n(^Xi) is the mean 
galaxy density at the distance Xi of the i'th galaxy. We determine n(^Xi) by integrating our estimate 
of the observed luminosity function, allowing for the effects of magnitude errors as described in 
Paper 1 (eq. [3]). Histograms of these predicted distributions are plotted as the dotted lines in 
Fig. 1. 

To apply this weighting scheme, we need a model for Js^s). As we showed in LEPM, large-scale 
clustering in the Stromlo-APM survey, and in other surveys, is well described by the linear power- 
spectrum of an Qh = 0.2 scale-invariant CDM model, hence we have calculated weights from (1) 
using this model. We find almost identical estimates of ^(s) if J3 for an Qh = 0.5 CDM model is 
used to calculate the weights. The estimates of ^(s) on large scales are not sensitive to the weighting 
scheme provided realistic values for J3 are used. 

In order to allow for the survey boundaries and selection function, we generate a random catalogue 
which fills the same volume as the galaxies and has the same selection function. Random points 
are generated within the distance range Xmin to x^ax according to a selection function obtained 
by integrating the observed luminosity function (Paper 1, eq. [11] and Fig. 4b). We choose 
Xmin = 5/i~-^Mpc, x^ax = 400/i~-^Mpc and a ratio of random points to galaxies Nj./Ng ~ 10. The 
sky coordinates of the points are chosen from a uniform distribution over each survey field. 



The standard estimator for ^(s) is 



Ug Wgr{s) 



where Wgg(^s) is the summed product of the weights (1) of galaxy pairs in separation bin s, Wgj.{s) 
is the equivalent quantity for galaxy-random pairs and Ug and Uj. are the mean densities of galaxies 
and random points, calculated with a minimum-variance estimator (Paper 1, eq. 4). The problem 
with this estimator is how one copes with fluctuations in the galaxy density Ug. By using a large 
enough random sample, fluctuations in Uj. can be made negligible, but for real redshift surveys, the 
actual density Ug for a subsample may not correspond to the expected density Ue for a homogeneous 
Universe. See, for example, Davis et al. (1988) and Maurogordato and Lachieze-Rey (1991) for a 
discussion of this problem. 

Hamilton (1993) has pointed out that one can measure ^(s) independently of any assumed galaxy 
density with the following estimator 

Here Wj.r{s) is the summed product of the weights of random-random pairs. Note that the relative 
densities of galaxy and random points measured at separation s are automatically accounted for 
by this estimator — there is no need to assume an overall galaxy density Ug. In this respect this 
estimator is similar to the 'ensemble' estimator of the angular correlation function w(^9) measured 
from counts of galaxies in cells 

These estimators are affected only to second order by density fluctuations related to the sample 
boundaries, whereas 'direct' estimators are affected to first order by such density fluctuations. 
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Of course, when comparing ^ between subsamples of a catalogue, one must be wary of the effects 
of such fluctuations in galaxy density because different parts of space are probed by the different 
samples. For example Fig. 1 shows that the overlap in volume between samples (b) and (d) is 
rather small, and that most of the galaxies are from independent volumes. Table 1 shows that the 
actual galaxy density Ug varies by up to 15% from the expected mean density Ue (given simply 
by integrating the luminosity function over the appropriate magnitude range) for the luminosity- 
selected samples, and so the observed behaviour of ^(s) determined with (2) at small amplitudes 
must be interpreted with some caution. One might choose to use the expected density Ue in equation 
(2) rather than the actual density Ug when normalising 1 + ^. However, this procedure leads to a 
positive 'tail' in ^ for over-dense samples and a negative 'tail' for under-dense samples, thus making 
comparison between samples difficult. By using the density-independent estimator (3), one does 
not have to assume a density for calculating ${s), normalisation is determined automatically from 
those galaxies at each given separation. Hence the density-independent estimates should be much 
more reliable than the density-dependent estimates. 

In Fig. 2 we plot ^(s) measured from the Stromlo-APM survey using both density-dependent (2; open 
symbols) and density-independent (3; solid symbols) estimators. In Fig. 2a we see that the density- 
dependent estimator finds considerably more power on large scales than the density-independent 
estimator, even when applied to the whole survey for which Ug is relatively well determined. The 
difference in estimates is due to the slight mis-match between the radial density functions of the 
galaxies and the random catalogue (Fig. la) — there is a slight underdensity in galaxies on scales 
^ 200/i~"^Mpc compared with what we would expect from the best-fit Schechter luminosity func- 
tion. This slight mis-match is detected by the density-dependent estimator as increased large-scale 
clustering. The density-independent estimator is much less sensitive to large-scale gradients in the 
data, in this case large-scale gradients in the relative galaxy/random density. 

A possible concern with the density-independent estimator is that it is removing intrinsic large-scale 
power in the galaxy distribution, not just artificial gradients due to uncertainty in the selection 
function. We can address this concern by analysing a volume-limited subsample of the catalogue, 
since one does not need to know the selection function to analyse a volume-limited (i.e. uniform 
density) sample. Additionally, no variable-weighting scheme is necessary for such a sample. Fig. 3 
shows the redshift-space two-point correlation function measured from a sample volume-limited to 
200/i"^Mpc (428 galaxies brighter than Mbj = -19.7). We now see very little difference between 
the density-dependent and density-independent estimators, suggesting that the density-independent 
estimator is not seriously 'filtering out' large-scale power. In fact, the density-independent estimator 
shows very slightly en/iancec? large-scale power over the density-dependent estimator for this sample. 

While the volume-limited sample shows slightly more large-scale power than the density-independent 
estimate from the full sample, it does not show the very large-scale clustering at s ~ 100/i~"^Mpc 
apparently detected by the density-dependent estimate from the full sample. Note the density- 
dependent estimate (Fig. 2a, open symbols) shows more large-scale power than our earlier determi- 
nation of ^(s) (LEPM, Fig. 3). For this earlier analysis, only galaxies within 300/i~-^Mpc distance 
were used, and the selection function was obtained by integrating the best-fit pure Schechter lu- 
minosity function. As discussed in Paper 1, a Schechter function convolved with a Gaussian helps 
correct for random magnitude errors in the data and provides a better fit to the observed luminosity 
function. It is the combined effect of a slightly different selection function and including galaxies 
beyond 300/i~-^Mpc which yields the increased large-scale power found with the density-dependent 
estimator in Fig. 2a. The density-independent estimator in Fig. 2a is in good agreement with the 
earlier (density-dependent) estimate in LEPM (which barely changes if a density-independent esti- 
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mator is used). The lack of sensitivity to the limiting distance and exact form of selection function 
is an important advantage of the density-independent estimator over the density-dependent one. 

3.2 Comparisons Between Galaxy Samples 

The clustering measurements using both density-dependent and density-independent estimators for 
all six galaxy samples listed in Table 1 are presented in Figure 2. The dotted lines in Fig. 2 show 
^(s) predicted by F = 0.2 biased CDM linear theory to aid comparing the samples. 

The density-dependent estimates show large variations on scales s > 20/i~-^Mpc, with a trend of 
increasing large-scale power with luminosity. As discussed in the previous section, the density- 
dependent estimates are sensitive to errors in the estimated galaxy density Ug and we believe this 
apparent strong trend is caused mainly by changes in the sample volume, and is not an intrinsic 
luminosity dependence effect. 

The density-independent estimates show much smaller differences in galaxy clustering with luminos- 
ity; the low-luminosity subsample (b) is slightly less strongly clustered than the other subsamples; 
there is no significant difference between the middle (c) and high (d) luminosity samples. Though 
these variations are smaller than those seen between the density-dependent estimates, the stabil- 
ity of the density-independent estimates means that they are much more significant. The slightly 
steeper and lower amplitude correlation function of faint galaxies compared to bright galaxies seems 
to be a real luminosity dependence effect. 

In Figure 2 we also present ^(s) measured from samples (e) and (f), early and late type galaxies. 
As expected, the early types show significantly stronger clustering than late types. Due to a bias 
against classifying galaxies at large distances in the survey as early-type (reflected by the low 
value of {V/Vmax) in Table 1 of Paper 1), rather than generate the random N(^x) distribution 
from the measured luminosity function for galaxies of the appropriate type, we have instead fitted 
a fourth-order polynomial to the observed N(^x) (Fig. 1). For this reason, no estimates of the 
observed/expected density ratio Ug/ue are given in Table 1 for these two subsamples. The 'tail' in 
the density-dependent estimate of ^(s) ~ 0.1 for early type galaxies is almost certainly due to the 
low {V/Vmax) for this sample. 

Power-law fits, ^(s) = (s/so)~'^% to the density-independent estimates over the range 1.5-30 /i~-^Mpc 
are given in Table 1 for each sample. For sample (a), all galaxies, we find a power-law index 
7s = 1.47, shallower than 7 ~ 1.7 measured in real space {e.g. Davis and Peebles 1983, Bean et 
al. 1983) and as determined from the angular correlation function w{9) {e.g. Groth and Peebles 
1977, Maddox et al. 1990c). This difference is due to redshift-space distortions (cf. the following 
section). Our estimates of 7^ and sq are slightly less certain than earlier determinations since our 
sparse-sampling strategy was designed to minimise errors in ^ on scales s ~ 20/i~"^Mpc, where the 
amplitude of ^ is low, not on small scales s ~ sq. 

The power-law index 7^ becomes progressively steeper for lower luminosity galaxies, changing from 
1.41 for the brightest sample to to 1.80 for the faintest sample. The correlation length for the 
lowest luminosity galaxies is sq = 4.9/i~-^Mpc, slightly lower than for higher luminosities which have 
So ~ 6/i~"^Mpc. This corresponds to a factor of 1.7 in the correlation amplitude at 10/i~"^Mpc, but 
redshift-space distortions mean that this cannot be interpreted directly in terms of the relative bias 
factors. 
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Early type galaxies show a very shallow 7^ ~ 1.25 due to the large amplitude of ^ on scales ~ 10- 
30/i~"^Mpc. As mentioned in Section 4, we believe this shallow slope is due to redshift distortions. 
The correlation length sq = 9.6 is significantly larger than for late type galaxies, corresponding to a 
ratio of 2.4 in amplitude at 10/i~"^Mpc; again redshift-space distortions must be taken into account 
in order to relate this value to relative bias factors. 



4 The Real Space Correlation Function ^(r) 

The clustering results presented in the preceding section are of course affected by redshift-space 
distortions to uniform Hubble flow. On small scales, random peculiar velocities will cause clustering 
to be underestimated, while on large scales coherent bulk-flows will cause the clustering amplitude 
to be overestimated (eg. Kaiser 1987). Moreover, we expect that early type galaxies, which are pref- 
erentially found in high density regions, will be more strongly affected by redshift-space distortions 
than late type galaxies found in the field. 

In order to measure clustering unaffected by redshift-space distortions, one must somehow integrate, 
or project, over the radial distance coordinate. There are several ways of doing this. 

One can measure the angular correlation function w(^9) and invert it using Limber's equation. The 
problem with this method is that w(^9) measured from a redshift survey is noisy compared to 
the w(^9) that can be measured from a much larger photometric catalogue without the redundant 
redshift information. Additionally, the inversion is sensitive to the selection function for the relevant 
galaxy type. Due to the difficulty of classifying 17th mag galaxies on Schmidt plates, the luminosity 
functions for early and late-type galaxies in this survey are subject to large errors. As described in 
§4.1 we circumvent these problems by measuring w(^9) for the fully-sampled APM Bright Galaxy 
Catalogue (Loveday 1989) and estimating the selection function S{z) by smoothing the observed 
iV(z) distribution for galaxies of appropriate type in the Stromlo-APM survey. 

Another approach, followed by Davis and Peebles (1983), is to calculate the full redshift space 
correlation function ^(ct, tt) as a function of the two components of separation parallel (tt) and 
perpendicular (ct) to the line of sight, perform the integral iH!(c7) = ^(^a,'K)d'K and then to invert 
the resulting projected correlation function iH!(c7) to obtain ${r). However, since ^(ct, tt) is now 
calculated on a 2d grid, it suffers from shot-noise due to the small number of galaxy-galaxy pairs in 
each (ct, tt) bin and so is very noisy. This is especially serious for our sparse-sampled redshift survey 
on small scales, and so this method is not used here. 

A third method, when one has a sparse-sampled redshift survey drawn from a larger parent cata- 
logue, is to calculate the projected cross-correlation between the redshift survey and its 2d parent 
survey. This projected correlation function is easily inverted to give ${r). This method gives the 
most stable and reliable estimates of ^(r), and we apply it to our survey in §4.2. 

4.1 Inversion of lo(^) 

In constructing the Stromlo-APM survey, Loveday (1989) inspected all APM galaxy candidates 
brighter than bj = 16.57 and assigned each galaxy a morphological type. There are 4439 early 
type and 8844 late type galaxies in the APM Bright Galaxy Catalogue (APMBGC) and so we can 
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calculate w(^9) for early and late type galaxies much more accurately from the APMBGC than from 
the 1:20 sparse-sampled redshift survey. As an additional bonus, the brighter mag limit of the 
APMBGC (16.57 vs 17.15) means that galaxy typing should be more reliable and complete (only 
164 out of 13447 unmerged APMBGC galaxies are not classified as early or late type). Of course, 
redshifts, and hence luminosities, are not known for the vast majority of APMBGC galaxies, and 
so this analysis cannot be applied to the luminosity-selected samples. 

We have estimated w(^9) using the estimator 

where Ngg, Ngj. and Nj-r are the number of galaxy-galaxy, galaxy-random and random-random pairs 
at angular separation 9 and Aw is a correction for the integral constraint. 



/ / w{9i2)dnidn2, (6) 

J J survey 



Aw 

I survey 



(Groth and Peebles 1977). The correction Aw is estimated in practice by calculating w{^9) without 
the correction, integrating w{^9) over all elements of solid angle dQ,i in the survey area to obtain Aw 
and recalculating w{^9) with the correction added. A stable solution is rapidly reached by iteration. 

Fig. 4 shows w{^9) for all, early and late-type galaxies in APMBGC. We have fitted a power law 
w{^9) = A9^~~' from 0.1 to 5° to these estimates, with results shown in Table 2. We see that early- 
type galaxies have a steeper power law slope than late-types, in agreement with Davis and Geller 
(1976) and Giovanelli et al. (1986) but contrasting with the redshift-space measurements (Table 1). 
The integral constraint corrections Aw shown in Table 2 make a negligible difference to power-law 
fits on scales smaller than 5° but they do give some idea of possible systematic errors in the w(^9) 
estimates on large scales. 

We have used these power law solutions in the relativistic version of Limber's equation (Groth 
and Peebles 1977, Phillipps et al. 1978) assuming qo = 0.5. The selection function S{z) used 
in Limber's equation was determined separately for each galaxy type by smoothing the observed 
iV(z) for galaxies in the Stromlo-APM survey of the appropriate type and with bj < 16.57 with a 
gaussian of FWHM = 0.01. The resulting parameters tq and B for the spatial correlation function 
^(r) = {r/ro)~~' = Br~~' are shown in Table 2. We see that at r = l/i~-^Mpc, the clustering 
amplitude of early-type galaxies is more than a factor of three higher than that of late-type galaxies. 



4.2 Projected Cross-Correlation ^{(t) 
4.2.1 Method 

Probably the most reliable way of determining the real-space correlation function when one has a 
sparsely-sampled redshift survey drawn from a fully sampled parent catalogue is to calculate the 
projected cross-correlation function iH!(c7) between the redshift survey and its 2d parent catalogue, 

^{^) = r°° a\l^y^T^')dAy, (7) 

J — oo 

where the integral extends over all line-of-sight separations Ay for pairs of galaxies with constant 
projected separation a = y9 (^9 is the angular separation and y is the distance to the galaxy 
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in the redshift survey). This projected function can be directly inverted numerically to give a 
stable estimate of ${r), which is unaffected by redshift-space distortions. This method was used by 
Saunders et al. (1992) to measure ^(r) for IRAS galaxies using the QDOT redshift survey (Lawrence 
et al., 1994) and was earlier used by Lilje and Efstathiou (1987) to measure the cross-correlation of 
Lick galaxies with Abell clusters. 

Here, we cross-correlate the redshift survey samples listed in Table 1 with 36,276 galaxies brighter 
than bj = 17.15 in the APM Galaxy Survey. The resulting estimate of ^ uses all galaxies from the 
parent sample as well as the sparse subsample with measured redshifts and so minimises random 
errors and enables us to estimate ^(r) on much smaller scales than possible using the redshift survey 
galaxies alone. A further advantage of this procedure for our analyses is that the estimates of ^ are 
independent of uncertainties in the selection function for galaxies of a specific type. 

To estimate S we consider each redshift survey galaxy at known distance y, and count the number 
of APM galaxies Ngi^a) at projected separation a = yO and then compare this with the number of 
random points Nj.{(t) (scaled by the relative numbers of galaxies and random points) at the same 
separation. 



An estimate of the projected correlation function. 



AX„ (8) 



is thus calculated for each redshift survey galaxy, where the factor 

1 f°° 

p{(T,y)=^r—— ip{x)x^^{r)dx, = x^ + y^ - 2xy cos{a / y) (9) 
jy :i[a) Jo 

corrects for projection effects and biases introduced by assuming that ^(r) is negligible on scales 
r ^ y, ip{x) is the galaxy density at distance x and N is the surface density of galaxies in the 2d 
catalogue. The term AX^ is a correction for the integral constraint affecting each redshift survey 
galaxy. It is estimated by assuming a truncated power-law model iH!(c7) = (ct/ctq)"" for a < cimax-, 
zero otherwise, and integrating over all solid angle elements in the survey area. 



/ / S((Ti2)c?nic?n2. (10) 

J J survey 



AX, , 

/ survey 



We assume parameters a = 0.81, ctq = 165/i~"^Mpc and (Tmax = 20/i~"^Mpc, which give a reasonable 
fit to the final iH!(c7) obtained from the all-galaxy sample. The correction AXi varies with the 
distance y of the redshift survey galaxy, from ~ 0.1 at y ~ 300/i~"^Mpc to ~ 20 at y = 10/i~"^Mpc. 

Saunders et al. (1992) made an estimate Xi(^a) of iH!(c7) for each redshift survey galaxy at distance 
yi, and then formed a weighted average of the Xi to obtain a final iH!(c7). In order to mimimise the 
shot-noise in each Xi(^a) estimate, we instead chose to bin the redshift survey galaxies into distance 
bins of width Ay = 10/i~"^Mpc and then estimated Xi for each bin centred on distance yi. 

Estimating the Xi in an unbiased way requires the correct values for p(c7, y) and AXi (Eqn. 8), which 
in turn require the prior knowledge of ^(r) and iH!(c7) (Eqns. 9 and 10). Fortunately p and AXi are 
only weakly dependent on ^ and S and so an unbiased solution can be calculated by iteration from 
approximate initial estimates of ^ and S. We assume initial power-law forms for ^(r) and iH!(c7) and 
calculate Xi for each distance bin yi. Then we take a weighted average of the Xi's using weights 
designed to give the minimum-variance estimate of S (Saunders et al eq. 16). This estimate of iH!(c7) 
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is then numerically inverted using Eq. 26 of Saunders et al to obtain an estimate of ${r). Note 
that this inversion does not assume a power-law form for ${r), an assumption which can lead to a 
systematic bias in the slope 7 (Saunders 1994). Power laws are separately fitted to the projected 
and spatial correlation functions over scales 0.2-10 /i~-^Mpc and the process is iterated to obtain a 
stable solution. Typically five iterations are required. 

4.2.2 Test of the Method 

In order to test the method, we have used one of the CDM-like N-body simulations of Croft and 
Efstathiou (1994). These simulations combine a large volume (box length = 300/i~-^Mpc) with a 
spatial resolution of ~ 80/i~"^kpc, and so can be used to generate reasonable approximations to our 
redshift survey. We analysed a simulation from ensemble B, a model with non-zero cosmological 
constant A = A/[3Hq) = (1 — fig) = 0.8. Weights were assigned to each particle using the peak 
background split algorithm (Bardeen et al. 1986, White et al. 1987) and 'galaxies' selected within 
the APM area and with the Stromlo-APM selection function (Paper 1). 

In Fig. 5 (solid symbols), we plot the correlation function ^(r) measured directly from the real-space 
positions of the 34,120 'galaxies' in the simulation. Over separations 0.2-20 /i~-^Mpc, ^(r) is very 
well fit by a power-law ^(r) = {r/ro)~~' with 7 = 1.89 and tq = 6.2/i~-^Mpc. 

We then generated five sub-samples from this simulation. For each sub-sample, galaxies were 
selected from the full simulation at random with probability 0.05, the same sparse-sampling rate used 
for the Stromlo-APM Survey. Each of these five mock redshift surveys was cross-correlated, using 
the redshift-distance information, with the sky-projected data from the full simulation. Estimates 
of iH!(c7) and thence ^(r) were determined for each mock redshift survey using the method described 
above. In Fig. 5, the average of the ^(r) estimates from the cross-correlation procedure is shown 
by the open symbols and the error bars show the scatter between estimates. It can be seen that 
this estimate of ^(r) gives excellent agreement with the 'true' ^(r) from the full simulation. The 
increased noise on scales r > 20/i~-^Mpc is not surprising since the power-law fits to iH!(c7) and ^(r) 
used in (9) were only made over the range 0.2-10 /i~-^Mpc. A final power-law fit to the average ^(r) 
gives 7 = 1.94 ± 0.08 and tq = 6.1 ± 0.4, consistent with the direct determination of ^(r) from the 
full simulation. 

4.2.3 Comparisons Between Galaxy Samples 

Our estimates of ^(r) obtained by inversion of the projected cross-correlation function between all 
APM galaxies in the parent sample and each of the subsamples of the redshift survey listed in 
Table 1 are shown in Fig. 6. Power-law fits over the range 0.2 < r < 20/i~-^Mpc are shown by 
dashed lines and the values of the slope 7^. and correlation length tq are listed in Table 1, where 
again the quoted errors are obtained from the variance between nine bootstrap resamplings of the 
redshift survey sample. To aid comparing the samples, the dotted line in each panel of Fig. 6 shows 
^(r) predicted in a P = 0.2 biased CDM model. 

For sample (a), all galaxies, we determine a correlation function slope 7^. = 1.71 ± 0.05 and a 
correlation length tq = 5.1 ± 0.2, in good agreement with earlier determinations {e.g. Davis and 
Peebles 1983). 
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We see no significant difference between the clustering of X* galaxies (sample c) and brighter galaxies 
(sample d). This is perhaps a surprising result given some previous work {e.g. Hamilton 1988), but 
is consistent with our density-independent estimates of ^ in redshift space (Section 3). We find that 
intrinsically faint galaxies (sample b) have a much steeper correlation function slope, jj. = 2.09, and 
smaller correlation length, tq = 3.2, compared to L* and brighter galaxies (samples c & d). This 
comfirms the trends which are seen at lower significance in our redshift-space estimates of Section 3. 
The correlation amplitude for sample (b) is a factor ~ 2 lower than (c) and (d) at l/i~-^Mpc and the 
factor increases to ~ 4 at 10/i~-^Mpc. A steeper slope might be expected if low-luminosity galaxies 
are found mostly in cluster environments and thus co-habit with early-type galaxies, as suggested 
by the steeper luminosity function faint-end slope found in cluster environments compared with the 
field (eg. Binggeli et al. 1988). 

We see that for early-type galaxies (sample e) the slope is jj. = 1.85, steeper than for late-type 
galaxies (sample f) which have jj. = 1.64. These slopes and also the corresponding scale lengths are 
in good agreement with observations of w{9) (Table 2). This confirms that the very shallow slope 
(7s = 1.25) seen in the redshift space estimate of ^ for early type galaxies is due to redshift-space 
distortions. 

It is important to remember that the ^(r) estimates shown in Fig. 6 are for the cross- correlation of 
galaxies of specified type with galaxies of all types. The differences between the cross-correlation 
functions of the samples should be smaller than the differences between the auto- correlation func- 
tions. For example, suppose that the mean galaxy bias relative to mass is b, and that early and late 
type galaxies have bias be and 6/. The cross-correlation functions ^eg and ^ig will have amplitudes 
Beg oc be X b and Big oc bi X b, with ratio be/bi. The auto-correlation functions ^ee and will have 
amplitudes Bee oc b^ and Bn oc bf, with ratio {bejbiY. From Table 1, the ratio ieglilg ~ 2.3 at 
1/i-^Mpc. From Table 2, we see that the ratio ieelill ~ 3.5 and not 5.5 as expected. Formally 
this discrepancy is marginally significant compared to the estimated errors, but note that the mag- 
nitude limit used to estimate w is different to that used for S and so different volumes are being 
sampled. Also the amplitude B is sensitive to the power-law slope 7 used in the Limber inversion. 
As discussed by Saunders (1994), the fact that w{9) is not a pure power-law, but contains a break, 
can result in a systematic bias in the 'power-law' slope 7. 

4.2.4 Effect of Sampling Fluctuations 

One might argue that the observed segregation at low luminosities and lack of segregation at high 
luminosities could be due to sampling fluctuations, since, however ingenious our estimator, low- 
luminosity galaxies are necessarily closer to us than more luminous galaxies. Ideally, one would like 
to form a volume-limited sample containing galaxies over a range of luminosities at uniform radial 
density, and so calculate the clustering of galaxies of different luminosities within the same volume. 
Unfortunately, true volume-limited samples drawn from the Stromlo-APM survey are very small 
indeed. This is because of the difficulty of photographic calibration brighter than 15th magnitude 
— a strictly volume limited sample must have a minimum distance and upper luminosity limit, as 
well as a maximum distance and lower luminosity limit. 

By way of compromise, we have calculated the spatial cross-correlation functions for the luminosity- 
selected samples with additional distance limits imposed. In Fig. 7a, we plot ^(r) for the faint (solid 
symbols) and middle (open symbols) luminosity galaxies, limiting to galaxies in the range 70-140 
/i-^Mpc (the range of overlap for these two samples — see Fig. 1). We see that even when drawn 
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from the same volume, middle-luminosity galaxies still show stronger clustering than low-luminosity 
galaxies. In Fig. 7b, we compare the clustering of middle (solid symbols) and high (open symbols) 
luminosity galaxies in their range of overlap (90-220 /i~-^Mpc). No clear difference in clustering is 
seen, in agreement with Fig. 6, in which no explicit volume constraint was imposed. Finally, in 
Fig. 7c, we plot the clustering of middle-luminosity galaxies in the two different volumes: 70-140 
/i~-^Mpc (solid symbols) and 90-220 /i~-^Mpc (open symbols). It is indeed encouraging that the 
same luminosity galaxies have the same measured clustering in two different (albeit not entirely 
independent) volumes. We thus believe that the observed segregation at low-luminosities and lack 
of segregation at high-luminosities is a genuine effect, and is not due to sampling fluctuations in our 
survey. 

4.2.5 Separating Luminosity and Morphological Segregation 

One might also ask whether the differences we see between the clustering of faint and middle 
luminosity samples are in fact due to true luminosity segregation or just a different balance of 
morphological types at different luminosity. Conversely, if, as Hasegawa and Umemura (1993) 
claim, early and late type galaxies show luminosity segregation of opposite sign, then the lack of 
segregation between the middle and bright luminosity samples could be due to cancelling of effects 
for early and late types. 

Since our ^(r) cross-correlation estimates in Fig. 6 have such small error bars, it is worth investigating 
the dependence of clustering on morphology and luminosity separately, ie. by further dividing the 
early and late type galaxies by luminosity. The cross-correlation ^(r) measured for these new 
samples are shown in Fig. 8 and the results of power-law fits from 0.2 to 20/i~"^Mpc are given in 
Table 3. 

Both early and late type galaxies separately show evidence of luminosity segregation between the 
faint and middle luminosity samples. The signal for the faintest ( — 19 < M < —15) early-type 
galaxies goes negative on scales K 6-20/2,-^ Mpc and so for this subsample, the power-law fit was 
truncated at 5/i~"^Mpc. Evidence for luminosity segregation brighter than L* is marginal at best, 
although possibly the brightest late type galaxies show slightly enhanced large-scale clustering over 
L* galaxies. The steeper, higher amplitude clustering of early-type compared to late-type galaxies 
occurs for all luminosity classes except for the faintest one, which gives a rather noisy correlation 
function for early types. 

The results presented in fig. 8 thus show that luminosity and morphological segregation are both 
real, independent effects. 

5 Discussion and Conclusions 

We have presented estimates for the correlation function of various galaxy samples in redshift space 
using two different estimators. The differences in redshift space correlation functions determined 
using density-dependent and density-independent estimators highlight the problems in trying to 
determine clustering from a sample in which the mean density is not well-defined. We have shown 
that the density-independent estimator (3) provides a more reliable determination of galaxy clus- 
tering when analysing subsamples of a catalogue in which the actual galaxy density differs from the 
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expected density, or when the selection function determined from the observed luminosity function 
does not provide a perfect fit to the observed radial density. 

We have seen that the redshift-space correlation function is significantly affected by peculiar veloc- 
ities and present estimates of ^(r) unaffected by these distortions. We find that early-type galaxies 
show a steeper correlation function slope and larger correlation length than late type galaxies. Low- 
luminosity galaxies exhibit a steeper slope and smaller correlation length than L* galaxies, but no 
significant difference is seen in the clustering of L* and super-X* galaxies. 

Our results concerning morphological segregation of galaxies are consistent with earlier investiga- 
tions by for example Davis and Geller (1976), Giovanelli et al. (1986) and lovino et al. (1993), 
who all find that early type galaxies are significantly more strongly clustered, and with a steeper 
correlation function slope, than late type galaxies. 

Given our results for variation of clustering strength with luminosity, it is not too surprising that 
previous analyses using smaller samples have not all agreed on the existence of luminosity segrega- 
tion. Our results are consistent with the majority of analyses which found no luminosity segregation 
or only a small difference in the clustering of faint and bright galaxies (see references in §1). 

Our results are not consistent with those of Hamilton (1988) who finds significantly enhanced clus- 
tering of the brightest galaxies compared with L* galaxies. Hamilton devised a test for luminosity 
segregation insensitive to variations in galaxy density by comparing clustering of galaxies of different 
luminosity in the same volume. The correlation function as a function of absolute magnitude is built 
up by multiplying ratios of correlation functions measured in successive volume-limited samples of 
the data. Unfortunately, this technique also accumulates errors in ^ as one works away from the 
fiducial luminosity, and so it is hard to assess the significance of the apparently enhanced cluster- 
ing of the most luminous galaxies seen by Hamilton. Interestingly, when Hasegawa and Umemura 
(1993) repeated Hamilton's analysis after correcting the CfA magnitudes for internal and galactic 
obscuration, the luminosity effect is much weakened. 

In most previous studies, low luminosity galaxy samples have been dominated by Virgo and the 
local supercluster, and so even if one allows for variation in galaxy density between samples (eg. 
Davis et al. 1988, Maurogordato and Lachieze-Rey 1991), one is still comparing clustering of faint 
galaxies in one small volume of the Universe with bright galaxies drawn from a much larger volume. 
Moreover, one expects galaxy peculiar velocities to be larger in high density regions such as the local 
supercluster, and hence redshift-space distortions may have had a stronger effect on low-luminosity 
galaxy samples than on high-luminosity ones. 

The significantly fainter apparent magnitude limit of the Stromlo-APM survey compared with earlier 
surveys means that our low-luminosity sample is drawn from a much larger volume of the Universe 
than was possible before — even our lowest luminosity sample has a median depth of ~ 100/i~"^Mpc. 
Therefore the statistical fluctuations on our clustering measurements for faint galaxies should be 
small, enabling reliable comparison with high-luminosity galaxy samples. Indeed, we have demon- 
strated this by comparing clustering of different luminosity galaxies in the same volume. 

The observed variation of galaxy clustering with morphological type and the observed weaker clus- 
tering of low-luminosity galaxies is what one would expect in biased galaxy formation scenarios. 
However, the lack of luminosity segregation at brighter luminosities is not compatible with some 
simple theories of biased galaxy formation, such as "natural bias" (White et al. 1987, Valls-Gabaud 
et al. 1989) in which galaxies preferentially form in the peaks of the underlying mass fluctuations. 
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The mean absolute magnitudes for our three luminosity subsamples are —18.4, —19.6 and —20.6. 
According to the White et al. model, and using the relation (1) between circular velocity and abso- 
lute blue magnitude given by White, TuUy and Davis (1988), one would expect enhancements in the 
amplitude of ^(r) by roughly a factor of 1.5 from each luminosity subsample to the next. Instead, 
we see an enhancement by factor 2-4 between the first two samples, and no significant difference 
between the second two samples. Conceivably, a closely related biasing model, but modified by a 
feedback mechanism at high luminosities, might be able to explain our observations. 

In order to compare the correlation functions for the various samples in real and redshift space 
more directly, in Fig. 9 we plot the real-space {${r), solid symbols) and redshift-space {${s), open 
symbols) correlation functions on the same plot. We have also re-binned the galaxy pair counts 
into coarser separation bins in order to reduce the error bars. For the all galaxies sample (a), we 
see a surprisingly small difference between ^(r) and ${s). This suggests a relatively low value for 
the quantity /3 = f2°'®/6 ~ 0.3, where Q is the cosmological density parameter and b is the bias 
parameter for optically selected galaxies. Further discussion of this topic is postponed until the 
next paper in this series (Loveday et al. 1995). 

Note that the ratio ^(•s)/^(r) in the linear regime for different samples will differ due to both 1) 
varying amplitude of redshift-space distortion with changes in the bias parameter bt for the different 
samples and 2) differences between the cross and auto-correlation functions. As discussed in §4.2.3, 
an auto-correlation function will scale as b^ whereas a cross-correlation function will only scale as 
6f. These two effects pull the ^(s) auto-correlation function in opposite directions relative to the 
^(r) cross-correlation function and so the large differences seen for low-luminosity galaxies (b) and 
early-type galaxies (e) are quite surprising. A quantitative analysis of this problem is again deferred 
to Loveday et al. (1995). 

It is intriguing that morphological segregation is strongest on scales < l/i~-^Mpc whereas luminosity 
segregation is strongest on scales > l/i~-^Mpc, hinting that the weaker clustering of low-luminosity 
objects may be a purely primordial (biasing) effect, but that morphological segregation may be 
enhanced by environmental effects, such as galaxy interactions and merging. Well motivated and well 
specified models of biasing as well as more observational data are needed to make further progress 
in understanding morphological and luminosity segregation, and thus providing an important key 
to unlocking the secrets of galaxy formation. 

Acknowledgements We thank Andrew Hamilton and Will Saunders for useful discussions. 
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Appendix: Optimal Weighting of Galaxies in a Redshift Survey 
for Estimating ( on Large Scales 



Imagine a homogeneous catalogue of galaxies with a well defined mean density n. Assuming that 
^(r) <C 1, then a realistic estimate of the error in ^ is given by 

1 + 4'KnJ3 

^^{r) ~ ^ . (11) 

where, 

J3(r)= r x^^{x)dx, (12) 
Jo 

and Np is the number of galaxy pairs used to estimate ^(r) for this separation bin (Peebles 1973, 
Kaiser 1986). Essentially, this model assumes that galaxies occurr in clusters of iVc = 1 + ^T^nJ^ 
members, and so the number of independent pairs in an estimate of ^(r) is given by Np/N^. 

Eqn. 11 will give an accurate estimate of 6^ for a homogeneous sample — e.g. an N-body simulation, 
but real redshift surveys are not homogeneous — the observed number density n of galaxies decreases 
with distance x from the observer. How should Eqn. 11 be applied in this case, and how should we 
weight each galaxy to minimise 6^ on large scales? 

We can estimate an optimum weighting scheme by considering galaxies in concentric shells about the 
observer, whose width is much larger than the correlation length tq. Suppose there are M 'cc-shells' 
of width I ^ To and volume Vi centred on the observer. If the density of galaxies in shell i is Ui, we 
expect ViUi galaxies in this cc-shell. Around each galaxy, neglecting edge effects, we estimate ^(r) 
in a spherical shell ('r-shell') of volume AV = 47rr^Ar. If r <C a;, then in the absence of clustering 
we would expect each r-shell to contain AVui galaxies. Therefore the expected number of pairs 
at separation r ± Ar/2 in cc-shell i is (Np-) = n^^VilSV . If we actually count iVp^ pairs, then our 
estimate of ^(r), which we denote ^e('')) is given by 

M 

\\Ur) = ^ , (13) 

i=l 

where Wi is a weight (to be determined) given to each galaxy in cc-shell i, the sums run over the M 
cc-shells, and we neglect pairs that cross shells. 

Now the error in ^ for shell i is 
and so from (11), 

SNp^ = {Np^fl\l + 47r7i,J3). (15) 

The variance in ^e('') is given by 



M 



Y.WtNp^l.l + ^'Kn^J^f 



2\ _ J=l 



M 



.i=l 



(16) 
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Substituting Np^ = njV^AV ( assuming ^(r) <C 1), and replacing the sums over cc-shells with inte 
grals, we get 

^ J W'^n{xY[l + 4'Kn{x)J3\^dV 



AV 



Jw^n{xy 



dV 



(17) 



where dV = ojx^dx for a redshift survey cone of solid angle oj, and the integration limits are the 



chosen distance limits for the redshift survey. Differentiating with respect to W, 



d{6e) 
dW 



if 



j W^n{xf[l + A'Kn{x)J^fdV J W^n{xfdV- 

j Wn{xfdV j W'^n{xf[l + A'Kn{x)J^fdV = 



and so 



W 



(18) 



(19) 



[1 + 47rTC(a;) J3] ' 

Of course, this optimal weighting requires prior knowledge of J^ir), but in practice W is only weakly 
dependent on J^ir) if it is large enough, and a stable solution may quickly be reached by iteration. 



Using this optimal weighting scheme, the estimated variance in ^ for small ^ is given by 



AV 



n[xfdV 
f 47rTC(a;) J3] 



(20) 
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Tables 



Table 1: Sample definitions and properties 





Type 








Ug/Ue 


Is 


•So 




''o 


a 


All 


-22 


-15 


1757 


0.991 


1.47 ± 0.12 


5.9 ± 0.3 


1.71 ± 0.05 


5.1 ± 0.2 


b 


Faint 


-19 


-15 


473 


0.965 


1.80 ± 0.18 


4.9 ± 0.6 


2.09 ± 0.13 


3.2 ± 0.6 


c 


Middle 


-20 


-19 


661 


1.154 


1.60 ± 0.22 


6.4 ± 0.8 


1.66 ± 0.06 


6.0 ± 0.4 


d 


Bright 


-22 


-20 


544 


0.952 


1.41 ± 0.31 


5.8 ± 1.6 


1.79 ± 0.11 


5.8 ± 0.3 


e 


E&SO 


-22 


-15 


336 




1.25 ± 0.33 


9.6 ± 1.2 


1.85 ± 0.13 


5.9 ± 0.7 


f 


Sp&Irr 


-22 


-15 


1062 




1.49 ± 0.21 


5.3 ± 0.4 


1.64 ± 0.05 


4.4 ± 0.1 



Notes. — 7s and sq are the power-law fit parameters to the correlation function in redshift space 
measured with the density-independent estimator over the range 1.5-30 /i~-^Mpc. jj. and tq are 
the real-space power-law parameters over 0.2-20 /i~-^Mpc determined from cross-correlation with 

the 2d APM survey (§4.2). 



Table 2: Angular correlation function results for all, early and late type galaxies in APMBGC 



Type 


7 


A 


Aw 


B 


''o 


All 


1.77 ± 0.03 


0.22 ± 0.01 


4.2 X 10-3 


19.0 ±0.8 


5.31 ± 0.15 


Early 


1.87 ± 0.07 


0.40 ± 0.03 


1.2 X 10-2 


46.0 ±5.0 


7.76 ±0.35 


Late 


1.72 ± 0.05 


0.18 ± 0.01 


3.4 X 10-3 


13.2 ±0.8 


4.49 ±0.13 



Note. — Power-law fits [w = A6^~~') were made over the range 0.1-5°. The integral constraint Aw 
is estimated from the observed w(^9). The amplitude, B, and corresponding scale length, tq are for 
the spatial correlation function inferred from inverting Limber's equation. 



Table 3: Results of power-law fits over 0.2-20 h -^Mpc for cross-correlation of joint morphology- 
luminosity selected samples with the 2d APM survey. 



Type 




Early Typ. 


e 




Late Type 




M^in Minax 






''o 






''o 


-19 -15 


78 


1.94 ± 0.22 


3.4 ± 0.8 


344 


2.01 ± 0.10 


2.9 ±0.4 


-20 -19 


102 


1.85 ±0.06 


6.4 ± 0.7 


403 


1.68 ±0.09 


4.9 ±0.2 


-22 -20 


153 


1.80 ±0.13 


7.4 ± 1.1 


312 


1.53 ±0.10 


6.2 ±0.4 
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Figure Captions 



Figure 1 N(^x) distributions for various galaxy samples as labelled. The dotted histograms show 
the predicted distributions from luminosity function estimates (a-d) and a fourth-order poly- 
nomial fit (e&f). 

Figure 2 The spatial correlation function ^(s) estimated from the Stromlo-APM redshift survey 
using density-dependent (open symbols) and density-independent (solid symbols) estimators 
for various samples as labelled. Error bars show the rms variance between nine bootstrap 
resamplings of the data. The dotted line shows the prediction of a T = 0.2 biased CDM model 
(LEPM), and is shown to aid in comparing the samples. 

Figure 3 The redshift-space correlation function ^(s) determined from a subsample of the survey 
data volume-limited to 200/i~-^Mpc. Open symbols show results from the density-dependent 
estimator, solid symbols results from the density-independent estimator. 

Figure 4 The angular correlation functions w(^9) for all (open circles), early-type (filled circles) 
and late-type (filled squares) galaxies in the APM Bright Galaxy Catalogue. The dotted lines 
show power-law fits from 0.1 to 5°. 

Figure 5 Test of the cross-correlation procedure. Solid symbols show ^(r) determined directly 
from a fully-sampled CDM-like N-body simulation; error bars are from the scatter between 
four zones. The open symbols show the average ^(r) determined from cross-correlating five 
1:20 random samplings of the simulation with the 2d information from the full simulation; 
error bars show scatter between the five estimates. 

Figure 6 The real-space cross-correlation function ^(r) determined by inversion from the projected 
cross-correlation function iH!(c7) for the galaxy samples listed in Table 1. Error bars show the 
scatter between nine bootstrap resamplings of the Stromlo-APM catalogue. 

Figure 7 Distance limited cross-correlation functions, (a) 70-140 /i~-^Mpc volume — low (solid 
symbols) and middle (open symbols) luminosity galaxies; (b) 90-220 /i~-^Mpc volume — mid- 
dle (solid symbols) and high (open symbols) luminosity galaxies; (c) middle-luminosity galax- 
ies in the two different volumes: 70-140 /i~-^Mpc (solid symbols) and 90-220 /i~-^Mpc (open 
symbols). 

Figure 8 The real-space cross-correlation function ^(r) determined by inversion from the projected 
cross-correlation function iH!(c7) for galaxy samples selected by both morphology and luminosity. 

Figure 9 Comparison of real-space clustering (solid symbols) and redshift-space clustering (open 
symbols). 
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